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Abstract. 

We present an edge-on Keplerian disk model to explain the main component of 
^\ • the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N . 

The brightness distribution and spectrum of the line of bright masers are successfully 
modeled with high amplification of background radio continuum emission along 
velocity coherent paths through a maser disk. The bend seen in the position-velocity 
diagram is a characteristic signature of differentially rotating disks. For a central 
mass of 3OM0, suggested by other observations, our model fixes the masing disk to 
have inner and outer radii of ~ 270 AU and ~ 750 AU. 
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o 

i-£h ■ 1. Introduction: disks in high- mass star formation regions 

1 ■ 

Disks are expected to form during protostellar collapse, and low-mass 

■ stars seem to provide good observational evidence for the existence of 
disks (e.g. Qi et al., 2003; Fuente et al., 2003). The situation is less clear 
for high-mass stars. While in several cases velocity gradients in massive 
star forming regions have been detected on large scales (> 10,000AU, 
e.g. Sandel et al., 2003), evidence for compact disks on scales < 1000 AU 
remains sparse. One possible example is IRAS20126-I-4104, studied by 
Cesaroni et al. (1999). 

Class II methanol maser emission, a signpost of high mass star for- 
mation (Minier et al., 2001), offers a potential indicator of disks since 
it often shows linear structures both in space and position - velocity 
diagrams (Norris et al., 1998; Minier et al., 2000). One of the most 
striking examples of a maser line in both space and velocity is found 
in NGC7538-IRS1N (Minier et al., 1998). We present here the first 
quantitative Keplerian disk analysis of this maser without invoking the 
assumption of a single radius. We find compelling evidence for the disk 
interpretation in this case. 
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Figure 1. Position - velocity diagrams of the 6.7 (top) and 12.2 GHz (bottom) 
maser data of NGC7538-IRS1 N . The outliers visible in the 6.7 GHz data were not 
considered for the fitting. 



2. Basic theory of the model and fit to the data 

Because of the remarkable agreement of the positions of the maser 
emissions at 6.7 and 12.2 GHz (2-3 mas), we draw the conclusion that 
the maser is the result of the amplification of a background continuum. 
Consider than an edge-on rotating disk at a distance D from the Sun 
and a point at radius p = r/D along a path with displacement 6. The 
rotation velocity V{p) and its line-of-sight component v obey v/8 = 
V/p = O, the angular velocity. In Keplerian rotation f2 oc p~ 3 / 2 with 
o = £l(po) = Dy/GM/R^ where p Q = R Q /D is the outer radius and M 
the central mass. Assume a Gaussian frequency profile for the maser 
absorption coefficient and denote the radial variation of its magnitude 
at line center by the normalized profile rj(p) (J rjdp = 1). Then 



t(6, v) = T / r){p) exp 



' 5 V Av D 



where f5 = y/1 — (0/p)' 2 and dp/ f3 is distance along the path. The width 
A^d is taken as constant and we set it to be the spectral line width 
of the maser emission at 12.2 GHz at 9 = 0. In our case Ai>d= 0.4 
kms -1 . This value could vary with p because of temperature variation 
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and maser saturation. Saturation typically requires r > 10-15 across 
the disk radius, and the broadening of the absorption profile is then 
proportional to length in excess of this threshold (Elitzur, 1992). The 
saturation issue is not treated in this study. 

By comparing the maser brightness temperature T& at = with 
T b of the background continuum (10000 - 15000 K, Campbell, 1984) 
we find that r = 18.32 and 15.99 for the 6.7 and 12.2 GHz masers 
respectively. We parameterize the function rj(p) as oc p( _p ) and set 
p = —0.5, being the value that gives the best fits. The adopted distance 
to the source D is 2.7 kpc. 

Fitting. The free parameters in the fit are the angular velocity fi m 
and the ratio between inner and outer radii, h. We define fi m = l/2(Qj+ 
fi ), i.e. the mean of angular velocities at the inner and outer radii, 
as the angular velocity measured in the data (see straight slope in 
Figure 3). Because of the natural scalefreeness of the model, we have 
to give an estimate for the scale of the disk, e.g. the outer radius. We 
choose a value which corresponds to a central mass of 30 Mq, relying on 
previous estimates in NGC7538-IRS1 N (e.g. Campbell, 1984). It is to 
notice however that the model can produce equivalently good fits when 
this estimate is not taken into account. The model shows a lower limit 
at disk sizes corresponding to central objects of subsolar masses. It still 
performs well when setting the central object to 100 M , although then 
the outliers at 6.7 GHz cannot be part of the disk. 

We fit the position - velocity diagram of the 12.2 and 6.7 GHz 
methanol maser data observed using VLBI (the latter without outliers, 
see Figure 1). The 12.2 GHz data fit is considered the reference fit 
because of the higher spatial resolution in the data. 

The kink. The most important dynamical feature of the model is 
surely the bend in angular velocity appearing at | ^ |~ 17 mas, visible 
only in the position - velocity data at 12.2 GHz. This is naturally 
reproduced by our model and it can be explained in the following way. 
When \6 \< 0k the coherence path participating to the maser extends 
over the entire disk, leading to the highest r and so the brightest maser 
emission. When increasing the displacement above 0k, the longest co- 
herence path will be centered at velocities closer to the outer edge of 
the disk, making the increase in angular velocity with displacement 
angle (or velocity gradient) less pronounced. This will produce a clear 
bend in the line tracing the maximum optical depth r at every value of 
0. The value for Ok gets smaller and the bend stronger the smaller Avn 
or the smaller h. The presence of the bend is a unique signature for a 
differentially rotating disk (e.g. a Keplerian disk as the one we model 
here), since a solid-body rotation would show a linear velocity gradient. 
We consider the independent observation of the bend in data and model 
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Figure 2. Data (left) and model (right) of the position - velocity diagram of the 
optical depth r of the main spectral feature of NGC7538 IRS1 N at 12.2 GHz. In 
the right panel the values for Q m and h as well as an estimate of the goodness of 
the fit (X 2 ) are listed. The blue slope in both panels represents f2 m . 



to be the most compelling evidence for a circumstellar rotating disk in 
NGC7538-IRS1N. 

Because of the high requirements in sensitivity and spatial resolu- 
tion, as well as on the smooth distribution of the masing material in 
the disk, this bend was not seen when modeling the H2O maser disk in 
NGC4258 (Watson and Wallin, 1994). It has been detected in thermal 
line observations and modeled with the help of multiple dynamical 
components (e.g. Belloche et al., 2002). Our modeling confirms that 
there is no need for multiple dynamical components when trying to fit 
the bend. 



3. Results 



The results of the fitting of the 12.2 GHz maser data is shown in the 
figures 2 and 3, J7 m = 0.055 and h = 0.36. These values correspond to 
a disk where the masing methanol is confined between 270 and 750 AU 
from the central protostar, which has a mass of 30 Mq. The model 
can reproduce with high accuracy both the brightness profile and the 
spectrum of the maser emission as well as the indications for a bend 
in the velocity gradient. The fitting of the 6.7 GHz data gave similar 
results, i.e. Vl m = 0.055 and h = 0.36. Following the line of maximum r 
in the position - velocity diagram in the adopted model we can consider 
the outliers at 6.7 GHz to be part of the disk. 
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Figure 3. Brightness profile (left) and spectrum (right) of the main spectral feature 
of the methanol maser of NGC7538 IRS1 N at 12.2 GHz. Filled blue lines are the 
data, dashed red lines the model. 



4. Final remarks 



The present study has produced the up to date most accurate model of a 
circumstellar disk traced by maser emission at small scales ( < 1000 AU) 
around a massive protostar. NGC7538-IRS1 N seems to be a special 
case, where the masing material appears to be distributed extremely 
smoothly in the disk and not to be particularly turbulent. 

The existence of the outliers in the 6.7 GHz data however indicates 
that the disk might have a more clumpy structure when approaching 
the outer edge: the gravitational potential being weaker it is possible 
to imagine the disk to get disrupted giving rise to local random en- 
hancements of the optical depth offering suitable conditions for maser 
amplification. In fact the outliers seem to be very close to the line of 
highest optical depth in an edge-on disk predicted by our model. 

A further issue is to decide whether the background continuum am- 
plified by the inverted methanol in the disk comes from the known 
UCHII region in NGC7538-IRS1 N and illuminates the entire disk or 
from the central object solely. The latter case seems to be in accordance 
with photoevaporation studies, where the radio continuum would be 
responsible for both the methanol injection into the gas phase and 
the detected radio continuum outflow (Hollenbach et al., 1994). At the 
moment it is not possible to discriminate both scenarios. 

Finally it is important to notice that the present study of maser 
emission does not differ very much from more frequent studies of (op- 
tically thin) thermal line emission. Once taken the logarithm of the 
maser brightness map the problem reduces to the study of the value of 
the optical depth r, which is exactly the same as in the case of thermal 
emission. The advantage of maser emission is the possibility to achieve 
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high spatial resolution and get information about the dynamics at very 
small scales (e.g. 1 AU at lkpc). 
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